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Abstract 


Results of buckling and nonlinear analyses of the Space Shuttle external tank 
superlightweight liquid-oxygen (LO 2 ) tank are presented. Modeling details and 
results are presented for two prelaunch loading conditions and for two full-scale 
structural tests that were conducted on the original external tank. The results illus- 
trate three distinctly different types of nonlinear response for thin-walled shells sub- 
jected to combined mechanical and thermal loads. The nonlinear response 
phenomena consist of bifurcation-type buckling, short-wavelength nonlinear bending, 
and nonlinear collapse associated with a limit point. For each case, the results show 
that accurate predictions of nonlinear behavior generally require a large-scale, high- 
fidelity finite-element model. Results are also presented that show that a fluid-filled 
launch-vehicle shell can be highly sensitive to initial geometric imperfections. In 
addition, results presented for two full-scale structural tests of the original standard- 
weight external tank suggest that the finite- element modeling approach used in the 
present study is sufficient for representing the nonlinear behavior of the superlight- 
weight LO 2 tank. 


Introduction 

The new era of cooperation in space between the 
United States and Russia has created a requirement for 
the Space Shuttle to reach a 51.6° high-inclination orbit. 
Currently, achieving this orbit requires the payload of the 
orbiter to be reduced by approximately 10000 lb. To 
recover most of the lost payload, National Aeronautics 
and Space Administration (NASA) is developing a new 
lightweight external fuel tank for the Space Shuttle made 
primarily from an aluminum-lithium alloy. This new 
design, referred to as the superlightweight external tank 
(SLWT), is expected to weigh approximately 58 000 lb, 
which is approximately 8000 lb lighter than the light- 
weight aluminum external tank currently in service. 

An important consideration in the design of the 
SLWT is the nonlinear behavior of the thin-walled 
regions of the structure that experience compressive 
stresses. Local or global buckling of the shell wall can 
cause the thermal protection system (TPS) to separate 
from the tank, which could cause the vehicle to fail. As 
part of a plan to ensure that the design does not have a 
shell-wall instability response, accurate predictions of 
the nonlinear response of the SLWT are needed. Accu- 
rate predictions of the nonlinear response of the SLWT 
require a large-scale, high-fidelity finite-element model 
to represent the complex structural details of the SLWT 
and a robust nonlinear shell analysis capability that can 
predict local and general instability buckling modes. 

One component of the SLWT that experiences sig- 
nificant compressive stresses is the liquid-oxygen (L0 2 ) 
tank. A primary goal of this paper is to present results 
that give an indication of how much additional load, 
beyond the operational loads for two critical prelaunch 
loading conditions, the L0 2 tank can withstand before 
buckling or exhibiting severe bending gradients. Toward 


that goal, this paper presents results of nonlinear analyses 
of the L0 2 tank conducted by NASA Langley Research 
Center personnel for the two prelaunch loading condi- 
tions. In addition, results arc presented for two full-scale 
structural tests of the original L0 2 standard-weight tank 
(SWT) tank configuration, which differs primarily in 
wall thickness, weight, and material type from the SLWT 
L0 2 tank. These results demonstrate that the finite- 
element modeling approach used in the present study 
yields accurate representations of the nonlinear structural 
response of the SLWT L0 2 tank. 

A second goal is to present results that illustrate 
some types of nonlinear behavior, and the sensitivity of 
those behaviors to initial geometric imperfections, that 
may be encountered in the design of liquid-fuel launch 
vehicles. Thus, results for the two prelaunch loading con- 
ditions and the full-scale structural tests are presented 
that illustrate three distinctly different types of nonlinear 
responses of thin-walled shells subjected to combined 
mechanical and thermal loads. The nonlinear response 
for the first prelaunch loading condition is characterized 
by a bifurcation-type buckling response that is insensi- 
tive to initial geometric imperfections and exhibits sub- 
stantial postbuckling load-carrying capacity. In contrast, 
the nonlinear response for the second prelaunch loading 
condition is characterized by a short- wavelength bending 
gradient that grows in amplitude in a stable manner with 
increasing load. The nonlinear responses for the two full- 
scale structural tests are characterized by a limit-point 
response of a doubly curved shell segment that buckles 
locally. The magnitude of the load corresponding to the 
limit point is shown to be sensitive to initial geometric 
imperfections in the L0 2 tank. 

For each of the three different nonlinear response 
phenomena presented, the modeling approach and imple- 
mentation used to develop a high-fidelity finite-element 



model that adequately represents the behavior are dis- 
cussed. An overview of the SLWT L0 2 tank and inter- 
tank structure and loading conditions is presented, the 
analysis code and finite-element modeling details are dis- 
cussed, and the results of the linear bifurcation buckling 
analyses and the nonlinear analyses are presented. 
Although the modeling approach and implementation 
presented are for a specific launch vehicle, they represent 
a general method that can be applied to the design of 
future liquid-fuel launch vehicles. 

SLWT Structure Overview 

The Space Shuttle consists of the orbiter, two solid 
rocket boosters (SRBs), and the external tank (ET), as 
shown in figure L The ET consists of a L0 2 tank, a 
liquid-hydrogen (LH 2 ) tank, and an intermediate struc- 
ture called the intertank (fig. 1). The intertank transmits 
the weight of the fuel, the ET structural weight, and the 
orbiter weight to the SRBs prior to launch and transmits 
thrust loads from the SRBs and the orbiter to the ET dur- 
ing ascent. The SLWT L0 2 tank is a thin-walled mono- 
coque shell made primarily of 2195 aluminum-lithium 
alloy . The L0 2 tank is approximately 49 ft long and has a 
maximum diameter of approximately 27.5 ft, as indicated 
in figure 2. The L0 2 tank consists of a forward ogive sec- 
tion made from 8 gore panels, an aft ogive section made 
from 12 gore panels, a cylindrical barrel section made 
from 4 barrel panels, and an aft elliptical dome section 
made from 12 gore panels. The coordinate systems used 
to locate the elements of the L0 2 tank and the intertank 
also are shown in figure 2. The coordinates {XT, Y, Z ) are 
typically referred to as the global coordinate system of 
the ET, and axial positions along the tank are indicated 
by the coordinate value of XT in units of inches. For 
example, the junction between the forward and aft ogives 
is indicated by writing XT - 536.74 in. Cylindrical coor- 
dinates are also used and are given by (r, 0, XT), where a 
positive value of 0 is measured from the positive Z-axis 
toward the positive F-axis, as shown in figure 2. 

The L0 2 tank also has a forward T-ring and an aft 
Y-ring frame that support a slosh baffle that prevents the 
fuel from sloshing during ascent. The slosh baffle, a 
lightweight (approximately 455 lb), thin-walled struc- 
ture, is supported by two deep, thin-walled rings at each 
end that attach to the forward T-ring and the aft Y-ring 
frame. Other parts of the L0 2 tank include a nonstruc- 
tural nose cone, a forged forward ogive fitting and cover 
plate, an aft spherical dome cap that contains the L0 2 
suction fitting and a covered manhole, and a vortex baf- 
fle attached to the base of the aft dome cap. The L0 2 tank 
gore and barrel panels are stretch formed, chemically 
milled, and then welded together. The panels are fabri- 
cated with substantial thickness tailoring to reduce struc- 
tural weight. The panels are made somewhat thicker at 


the welds to form a stiffener-like region that is used as a 
weld land. The primary role of the weld lands is to com- 
pensate for reduction in shell- wall strength that is caused 
by welding. Tapering the weld lands in thickness and 
width along their length reduces weight and alleviates 
stress concentrations in the shell that result from abrupt 
changes in thickness. 

The intertank is a right-circular cylinder that is made 
from 2090 and 7075 aluminum alloys and is shown in 
figure 3. The approximately 22.5-ft-long intertank has a 
diameter of approximately 27.5 ft and consists of six 45° 
curved panels that are stiffened longitudinally with exter- 
nal hat stiffeners and are referred to herein as skin- 
stringer panels. The intertank also has two massive 45° 
curved panels, referred to as thrust panels, located per- 
pendicular to the F-axis of the intertank (fig. 3) that are 
stiffened longitudinally with integrally machined exter- 
nal blade stiffeners. These eight panels are assembled 
into the intertank with mechanical fasteners and are 
attached to five large internal ring frames, a forward 
flange, and an aft flange, Longitudinal straps (referred to 
herein as roll ties) suppress lateral-torsional deflection of 
the ring frames. The main central ring frame, two thrust 
panel longerons, and the thrust panels are connected to 
each end of a tapered beam that is referred to herein as 
the SRB beam (fig. 3). The SRB beam spans the diame- 
ter of the intertank along the F-axis and has a maximum 
depth (in the XT direction) of approximately 43 in. at its 
midspan. Forged fittings (referred to herein as SRB 
thrust: fittings) that are incapable of transmitting 
moments are fastened to the ends of the SRB beam. The 
primary role of the thrust panels is to diffuse the large 
axial loads introduced by the SRB’s into the intertank 
and then into the LQ 2 tank shell wall. The SRB beam 
compensates for the eccentricity of the concentrated 
loads introduced by the SRB’s. The SRB beam also sup- 
ports the loads normal to the intertank (parallel to the 
SRB beam) at the SRB attachment points. The intertank 
also has a 46-in-high by 52-in-wide frame-reinforced 
nonstructural access door located along the cylinder gen- 
erator at approximately 0 = 146°. 

Critical Prelaunch Loading Conditions 

Several critical SLWT loading conditions have been 
identified by the members of the SWLT team at the 
NASA George C. Marshall Space Flight Center and at 
the Lockheed-Martin Manned Space Systems Division. 
The general characteristics of the critical loads are indi- 
cated in figure 4, These loads consist of the wind or aero- 
dynamic pressure loads, the structural weight or inertia, 
the pressure exerted on the shell wall by the L0 2 , the 
ullage pressure inside the tank, the interface forces 
exerted by each SRB (indicated by the vectors Rj and R 2 
in fig. 4), the interface forces between the intertank and 
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the LH 2 tank (indicated by the vectors F and M in fig. 4), 
and the thermally induced loads associated with the cryo- 
genic temperatures. 

Two critical loading conditions addressed herein 
correspond to prelanch fueling conditions that occur 
when the Space Shuttle is on the launch pad. Prior to 
launch, the LH 2 tank is filled with LH 2 , and then the L0 2 
tank is filled with LO 2 . The first loading condition con- 
sidered corresponds to a full LH 2 tank and an empty L0 2 
tank. For this condition, there is no pressure and no tem- 
perature change in the L0 2 tank. However, the lower 
45 in. of the intertank are subjected to an ax i symmetric 
uniform through-the-thickness temperature field that var- 
ies linearly from -423°F, where the intertank is attached 
to the LH 2 tank, to 50°F at the top of the LH 2 tank for- 
ward dome. The nominal ambient temperature of the 
LO 2 tank and the intertank prior to fueling is 50°F. The 
SRB interface forces are given by 

R, = - 224.09-21 - 52.223 j - 28.954k kips 

R 2 = - 343,6241 + 48.261 j - 30.754k kips 

where the vectors Rj and R 2 are shown in figure 4 and 
where i, j, and k are standard orthonormal base vectors 
associated with the XT, Y, and Z axes, respectively. Simi- 
larly, the interface force and moment between the inter- 
tank and the LH 2 tank are given by 

F = 541.5931 + 9.614j + 63.494k kips 

M = -310.5001 + 10 715. 745 j + 16 828.589k in-kips 

The L0 2 tank is also subjected to wind pressure that has 
a resultant force given by -5.652J - 3.786k kips. 

The second prelaunch loading condition considered 
corresponds to full LH 2 and L0 2 tanks, but without 
ullage pressure. For this condition, the L0 2 tank is filled 
to XT = 447 in. (fig. 2), which corresponds to a depth of 
approximately 43 ft. For this case, the hydrostatic pres- 
sure distribution in the L0 2 tank is essentially axisym- 
metric, and the temperature distribution in the L0 2 tank 
and the intertank is axisymmetric and uniform through 
the thickness, as shown in figure 5. The specific weight 
of the L0 2 used in the present study is 0.04123 lb/in 3 , 
and the maximum hydrostatic pressure is approximately 
21 psig at the bottom of the tank. The corresponding 
weight of the L0 2 is approximately 1348 kips. 

The temperature field for the L0 2 tank shown in 
figure 5 has a uniform value of -297°F between 
XT — 447 in. and XT = 852.8 in,, but varies from a value 
of-150 9 F at the nose to -297°F at the free surface of the 
L0 2 (XT = 447 in.). The intertank is subjected to an axi- 
symmetric uniform through-the-thickness temperature 
field that varies in a piecewise linear manner from 


-423°F to 50°F and then to -297°F as the intertank is tra- 
versed from XT- 1129 in. to XT - 852.8 in. The temper- 
ature field shown in figure 5 represents the changes in the 
temperature distribution that occur from the nominal pre- 
fueling temperature of 50°F. 

The SRB interface forces for the second prelaunch 
loading condition are given by 

Rj = -906.447i-105.605j- 30.502k kips 

R 2 = —1035.789i + 1 08.1 28 j - 32.557k kips 

Likewise, the interface force and moment between the 
intertank and the LH 2 tank are given by 

F - 568,113i + 3.014j + 66.913k kips 

M = -354.488i + 11 249.31 8 j + 19443,206k in-kips 

and the L0 2 tank is subjected to a wind load with a 
resultant force given by -5.537j - 3,854k kips. 

Analysis Code and Finite-Element Modeling 

The results of the linear bifurcation buckling and 
nonlinear analyses were obtained with the Structural 
Analysis of General Shells (STAGS) nonlinear structural 
analysis code for general shells (ref. 1). The finite- 
element models of the SLWT used in the present study 
are very complex and include many structural details and 
the skin thickness variations or tailoring used to reduce 
structural weight. STAGS was chosen for analyzing the 
SLWT because of its robust state-of-the-art nonlinear- 
equation solution algorithms and its general user-input 
capability that is convenient for modeling branched 
shells typically used for launch vehicles. In particular, 
STAGS uses both the full and modified Newton methods 
to obtain an accurate nonlinear solution, and large rota- 
tions in the shell are represented by a corotational algo- 
rithm at the element level. The Riks arc-length projection 
method is used to continue a solution past limit points. 
STAGS permits complex geometries, loading conditions, 
and initial geometric imperfections to be modeled in a 
direct manner by the use of user-written subroutines that 
are essentially independent of the mesh discretization. 
For example, these user-written subroutines allow the 
user to define reference surface geometries; tapered shell 
walls and stiffener cross sections; and complex nodal 
force, temperature, and pressure distribution functions in 
a direct manner, using the FORTRAN computer lan- 
guage. This feature greatly simplified the definition of 
the finite-element models and the mesh convergence 
studies conducted in the present study. A description of 
how the features of STAGS were used in the present 
study to model the SLWT L0 2 tank and intertank 
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appears in this section. Details of how the prelaunch 
loads are simulated are also presented. 

In modeling the SLWT L0 2 tank and the intertank, 
several assumptions were made to simplify the finite- 
element models. Limited parametric studies were made 
to determine the adequacy of the assumptions and simpli- 
fications, For each case, the modeling assumptions used 
to simplify the finite-element models were found to be 
acceptable for analyzing the nonlinear behavior of the 
SLWT L0 2 tank. 

L0 2 Tank Modeling Details 

One very useful method for finite-element-mesh 
generation in STAGS involves defining the number of 
rows and columns of nodes in a segment of a specific 
type of shell (reference surface) unit, such as a cylinder, a 
cone, or a sphere. STAGS includes a library of several 
different predefined shell units that are based on a spe- 
cific type of surface parametric representation. For 
instances where a different surface parametric represen- 
tation for one of the shell units in the STAGS library is 
needed, or when a type of shell unit is needed that is not 
in the STAGS library, the surface information can be 
input into STAGS with the user-written subroutine called 
LAME. 

For the forward and aft ogive sections of the SLWT 
LO 2 tank, a parametric representation defined in terms of 
the global circumferential coordinate 0 (defined in fig. 2) 
and a local axial coordinate of the shell unit were put into 
LAME. The ogive geometry is based on a nominal 
612.0-in. meridional radius of curvature and a 165.5-in. 
maximum polar radius of curvature. The reference sur- 
faces for the forward ogive fitting and cover plate indi- 
cated in figure 2 were input as an ogive segment and a 
circular segment of a plane, respectively. The ogive seg- 
ment has an axial (XT) length of 4 in. and polar radii 
of 25.676 in. and 28.950 in. at XT = 371 in. and 
XT = 375 in., respectively. Similarly, the flat reference 
surface of the cover plate has a radius of 25.676 in. 

The reference surface for the barrel section of 
the L0 2 tank was input as a right-circular cone with 
polar radii equal to 165.359 in. and 165.373 in. at 
XT= 748,67 in. and XT == 843 in., respectively, and the 
reference surface for the aft dome section was input as a 
truncated ellipsoid attached to a smaller spherical cap. To 
simplify the model, the L0 2 suction fitting and covered 
manhole were neglected in defining the reference surface 
(and mesh) of the spherical cap. This simplification is 
based on the reasoning that the local details of the spheri- 
cal cap will have a negligible effect on the behavior of 
the ogives and the barrel. The ellipsoid has a 1 65.50-in. 
semimajor axis in the radial direction and a 124.125-in. 
semiminor axis in the axial direction. The spherical dome 


cap has a 70.0-in. polar radius at XT - 951.526 in., a 
211.855-in. spherical radius, and a 11.9-in. axial length. 
The aft Y-ring section of the tank indicated in figure 2 
was modeled with four reference surfaces that connect 
the barrel to the aft dome and to the intertank. The refer- 
ence surface that connects the barrel to the intertank is a 
9.8-in-long, 165.373-in-radius cylinder. Two conical ref- 
erence surfaces connect the barrel to the aft dome: the 
first is an 8.0-in-long cone with outer and inner polar 
radii equal to 165.373 in. and 164.72 in. at AT = 843 in. 
and XT - 851 in., respectively, and the second is a 
4.273-in-long cone with outer and inner polar radii equal 
to 164.72 in. and 164.08 in. at XT —851 m. and XT = 
855.273 in., respectively. The fourth reference surface is 
a very short and essentially rigid cone that connects the 
aft end of the cylindrical surface to the intersection of the 
two conical sections. This cone was included in the 
model to represent adequately the stiffness of the Y-ring 
section where the shell segments join in order to prevent 
the Y-ring from passing through the adjacent shell walls 
during deformation. This modeling detail was also used 
to eliminate large artificial bending gradients in the 
Y-ring shell wall that caused convergence problems in 
the numerical solution of the nonlinear equations. 

The SLWT LQ 2 tank shell wall, chemically milled to 
reduce structural weight, has a highly variable thickness 
distribution. In the forward ogive, the thicknesses vary 
from 0.080 in. to 0.157 in. in both the meridional and cir- 
cumferential directions. Similarly, in the aft ogive and 
barrel sections the thicknesses vary from 0.081 in. to 
0.190 in. and from 0.140 in. to 0.385 in., respectively. In 
the aft dome the thicknesses vary from 0.088 in. to 
0.125 in. This complex thickness distribution and corre- 
sponding wall eccentricities were input into the STAGS 
user-written subroutine WALL with FORTRAN state- 
ments. In addition, the LO 2 tank is thickened locally in 
several regions that support an external cable tray, and 
these local thicknesses were included in subroutine 
WALL. 

Several meridional lines shown in figure 2 represent 
the weld lands of the gore and barrel panels that form the 
tank. These weld lands are substantially thicker than the 
adjacent shell wall, and their thicknesses vary in the 
meridional direction. The weld lands were modeled as 
discrete beams of rectangular cross section with variable 
width and depth. The variations in the weld-land cross- 
sectional geometry in the meridional direction were rep- 
resented by using a linear variation with the correspond- 
ing surface arc length. The variable properties of the 
beams were input into STAGS with the user- written 
subroutine CROSS. The circumferential location of 
each weld land was specified to be at the column of 
nodes in the shell unit that is closest to its actual location. 
This approximation becomes more accurate when the 
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circumferential mesh density is increased. Circumferen- 
tial weld lands are located where the forward and aft 
ogives are connected, where the T-ring attaches to the aft 
ogive and barrel, where the Y-ring attaches to the barrel 
and the aft dome, and where the elliptical part of the aft 
dome connects to the spherical cap. The weld lands 
located at the aft end of the barrel and elliptical dome 
were modeled as discrete rings for convenience, and all 
the other weld lands were modeled with shell elements. 
In addition, the two deep, thin-walled slosh baffle rings 
were modeled as discrete rings. The remainder of the 
slosh baffle was not included in the STAGS models 
because the stiffness and stability it provides to the slosh 
baffle rings are essentially included by modeling the 
slosh baffle rings as discrete rings. 

Parts of the forward ogive fitting and the entire cover 
plate contain integrally machined stiffeners. These struc- 
tures were modeled as an equivalent homogeneous ortho- 
tropic wall with the appropriate shell reference surface 
eccentricity. In addition, cutouts and other local details in 
the cover plate were neglected, and places in the ogive 
fitting and cover plate where the thickness is much larger 
than the nominal wall thickness were modeled using dis- 
crete rings. 

Intertank Modeling Details 

The reference surface for the intertank is a 
270.35-in-long, 165.373-in-radius right-circular cylinder. 
The six skin-stringer panels of the intertank consist of 
external hat stiffeners bonded to a variable-thickness 
wall. The skin thicknesses (including primary doublers) 
range from 0.067 in. to 0.221 in. The hat stiffeners range 
in thickness from 0.045 in. to 0.063 in., arc 2.50 in. deep, 
and taper in width from 2.57 in. at the shell wall to 
1.38 in. at the top of the hat. The hat stiffener spacing is 
7.20 in. The objective of the intertank model is to pro- 
vide an accurate representation of the overall intertank 
stiffness so that load can be transferred to the LQ 2 tank 
with as few finite-element degrees of freedom (DOF) as 
possible. To achieve this goal, the skin-stringer panels 
were modeled as an equivalent homogeneous orthotropic 
wall with a variable thickness that is eccentric to the shell 
reference surface. STAGS computed the equivalent 
homogeneous orthotropic wall properties by using the 
“smeared stiffener” capability of the user- written subrou- 
tines WALL and CROSS. This model includes the 
effects of all the hat stiffeners and the primary skin dou- 
bler plates. Eight extruded stringers fastened to four of 
the panels were not included in the model as discrete 
beams, but were approximated in a conservative manner 
as hat stiffeners with a maximum thickness of 0.63 in. 
The variable wall properties were input into STAGS with 
the user-written subroutine WALL. Cutouts in the pan- 
els, including the access door and its frame, were 


neglected based on the rationale that their influence on 
the transfer of load into the L0 2 tank is small. 

The two thrust panels of the intertank are blade- 
stiffened panels integrally machined from a single piece 
of aluminum alloy. A total of 65 skin thicknesses that 
range from 0.090 in, in areas away from the central 
region of the panels to 2.062 in. in the area next to the 
SRB thrust fittings was used to model the thrust panels. 
A wall thickness of 2.062 in. was placed in the cutouts 
that house the SRB thrust fittings to facilitate load diffu- 
sion from the SRB beam to the thrust panels. The two 
internal longerons fastened to the thrust panels and the 
SRB beam, however, were neglected. Each panel has 
26 variable- width blade stiffeners that are 2.06 in. deep. 
Thirty-six blade widths ranging from 0,180 in. to 
1.050 in. were used to model the thrust panels as an 
equivalent homogeneous orthotropic wall with variable 
thickness in a manner similar to the skin-stringer panels. 

The five internal ring frames and two (forward and 
aft) flanges of the intertank were modeled as discrete 
rings with variable and constant eross-secriona; proper- 
ties, respectively. For the ring frames, the section proper- 
ties were modeled in a piecewise manner with constant 
properties for each beam element that connects two adja- 
cent nodes. Because the discrete ring model of the ring 
frames does not account for lateral -torsional cross- 
sectional deformations, the effect of the roll ties shown in 
figure 3 is implicitly included. The SRB beam was also 
modeled as a discrete beam in which several beam ele- 
ments were used to simulate its variable depth in a piece- 
wise manner. The SRB beam, approximately 345 in. 
long, extends outboard of the shell reference surface to 
cause an eccentricity in the SRB interface forces of 
approximately 7 in., as shown in figure 4. Load diffusion 
from the ends of the SRB beam to the shell wall of the 
thrust panels is facilitated by connecting the node at each 
end of the SRB beam to the adjacent nodes on the thrust 
panels with high- stiffness beam elements. 

Load Simulation 

A primary goal of the SLWT L0 2 tank study pre- 
sented herein was to determine how much additional 
load, beyond the operational loads, the tank can with- 
stand before buckling or exhibiting severe bending gradi- 
ents that will damage the thermal protection system 
(TPS). The basic approach used in the present study to 
achieve this goal is to apply all of the loads illustrated in 
figure 4 to the model, except for the SRB interface loads. 
The nodes on the ends of the SRB beam, where the SRB 
forces act, were restrained so that the SRB interface 
forces became reactions and rigid body motion was elim- 
inated. Next, the applied loads were separated into two 
groups. The first group contained the LH 2 tank interface 
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force and moment, which were treated as the primary 
source of destabilizing compressive stresses in the L0 2 
tank that may occur at load levels greater than the corre- 
sponding operational load level. The second group of 
loads consisted of the L0 2 pressure (for the second load- 
ing condition), the wind load, the structural weight, the 
thermal load, and the weight of the slosh baffle located 
inside the barrel section of the L0 2 tank. Constant in 
value, the loads in the second group were part of the 
operational loads considered to be passive loads when 
determining the stability margin of safety of the L0 2 
tank. 

The simple L0 2 pressure distribution and the tem- 
perature distribution were input into the STAGS model 
with the user-written FORTRAN subroutines UPRESS 
and UTEMP, respectively. The substantially more com- 
plicated wind load required construction of a Fourier 
series representation of the pressure field that was input 
directly into subroutine UPRESS. The slosh baffle 
weight of approximately 455 lb was applied to the slosh 
baffle support rings at XT = 744.85 in. and XT = 851.0 in. 
(fig. 2) as eccentric, uniformly distributed line loads. The 
LH 2 tank interface force and moment were applied to the 
model with the least squares loading and moving plane 
boundary features of STAGS. The STAGS least squares 
loading feature used a least squares fit to convert concen- 
trated forces and moments applied at an axial location 
into statically equivalent shell-wall stress resultants. The 
moving plane boundary feature of STAGS enforced the 
geometric constraint that all nodes within the given plane 
remain coplanar during deformation. 

In performing linear bifurcation buckling and non- 
linear analyses with STAGS, two load factors, p a and p h , 
were assigned to the first (active) and second (passive) 
load groups, respectively. First, a linear analysis was 
conducted to verify that the SRB reactions calculated 
from the applied loads were reasonably close to the spec- 
ified values defined previously. For both prelaunch load- 
ing conditions, the i and k components of the reactions 
were in good agreement with the corresponding specified 
values. In contrast, the j components of the reactions 
were not in very good agreement with the specified val- 
ues. However, because the j components of the reactions 
act along the axis of the SRB beam, their effect is mostly 
contained within the intertank. As a result, the SRB beam 
reactions are reasonable approximations of the actual 
forces. 

For the linear bifurcation buckling analyses, the load 
factor for the passive load group was assigned a value of 
one (p b = 1), and the load factor for the active load group 
p a was defined as the eigenvalue. This approach pro- 
vided a linear prebuckling stress state in the model that 
was used in determining the eigenvalue. For the nonlin- 


ear analyses, the load factors for both load groups were 
increased simultaneously to a value of one (p a =p b = 1), 
which corresponds to the operational values of the loads 
and provides the proper nonlinear prebuckling state for 
the L0 2 tank. Then, the load factor p a of the primary 
destabilizing loads was increased until an instability was 
reached. 

SLWT Prelaunch Load Results 

Results are presented in this section for the two pre- 
launch loading conditions previously discussed. First, 
results are presented for the loading condition that has a 
full LH 2 tank and an empty L0 2 tank. Then, results are 
presented for the loading condition that has a full LH 2 
tank and a full L0 2 tank. For this second loading condi- 
tion, results presented were obtained from STAGS mod- 
els that neglect the slosh baffle ring stiffnesses. Other 
results obtained from models that include these ring 
stiffnesses show that neglecting the slosh baffle ring 
stiffnesses in the STAGS models yields similar structural 
deformations and conservative predictions of the L0 2 
tank load-carrying capacity at load levels greater than the 
level of the operational loads. 

Full LH 2 and Empty L0 2 Tanks 

Several different finite-element meshes were used in 
the study for analysis of the L0 2 tank subjected to the 
prelaunch loading condition with a full LH 2 tank and an 
empty L0 2 tank. As a first step toward identifying an 
adequate mesh with as few degrees of freedom as possi- 
ble, linear bifurcation buckling analyses were conducted. 
The passive loads associated with load factor p b were 
applied to the STAGS models as a linear prebuckling 
stress state (p b - 1), and the active (destabilizing) loads 
associated with load factor p a were used to obtain the 
minimum eigenvalue. 

The meshes investigated ranged from 104600 to 
213 500 DOF. The mesh identified as adequate for pre- 
dicting the linear bifurcation buckling behavior is shown 
in figure 6 and corresponds to 146700 DOF. This figure 
shows a buckle in the LO 2 barrel on the negative F-axis 
side of the tank. The reduction in mesh size from 213 500 
to 146700 DOF was done by increasing the mesh fine- 
ness in the local region containing the buckle shown in 
figure 6 and then by eliminating unneeded mesh refine- 
ment elsewhere, with care not to introduce spurious solu- 
tions. This step was facilitated by the use of the five-node 
and seven-node rectangular transition elements available 
in STAGS. The 104 600- and 146700-DOF models have 
the same general mesh arrangement shown in figure 6, 
but the level of local refinement of the mesh shown in 
figure 6 for the 1 46 700-DOF model is essentially twice 
that of the 104 600-DOF model. The eigenvalues for the 



104600' and 146700-DOF models are given by 
p a = 3.36 and p a - 3.26, respectively. These values corre- 
spond to loads that are approximately three times the 
magnitude of the operational loads. Based on the 
smoothness of the buckling mode shown in figure 6 and 
the 3-percent difference in the eigenvalues, the 
146700-DOF model is considered adequate for repre- 
senting the linear bifurcation behavior of the L0 2 tank 
for this loading condition; 

Next, the 104600- and 146 700-DOF models were 
used to conduct nonlinear analyses of a geometrically 
imperfect shell. An imperfection shape in the form of the 
corresponding linear bifurcation buckling mode and a 
negative value of the imperfection-amplitude-to-wall- 
thickness ratio A/t = 0.25 were used in these analyses. 
The thickness t in the ratio A/t, the minimum-gage wall 
thickness for the LO 2 barrel, has a value equal to 
0.140 in. In general, the sign of the linear bifurcation 
buckling mode is arbitrary and is determined by the spe- 
cific algorithm used to perform the calculations. For a 
shell structure with nonnegative Gaussian curvature, the 
sign of the eigenvector of the buckling mode indicates 
whether a specific region of the buckling mode of a 
curved surface is directed toward or away from the con- 
cave side of the surface; As a result, these two orienta- 
tions of the same buckling mode generally correspond to 
different degrees of nonlinear interaction between imper- 
fection shape and the membrane compressive stresses. In 
the study, the negative of the linear bifurcation buckling 
mode amplitude obtained from the STAGS models was 
used with the imperfection shape because it provided the 
strongest nonlinear interaction with the compressive 
stresses in the shell wall. The results of these analyses, 
shown in figure 7, are presented; as a plot of load factor 
p a versus the normal displacement of the shell wall at the 
largest crest of the buckle defined by the coordinates 
XT- 787.97 in. and 0 = 300.94°. (See fig. 2 for coordi- 
nate definitions.) The dashed and solid lines shown in the 
figure correspond to the 104 600-DOF and 146 700-DOF 
models, respectively. The unfilled square symbols and 
the filled circular symbols shown in the figure corre- 
spond to the actual analytical results obtained for the 
104600-DOF and 146 700-DOF models, respectively. 
The results shown in figure 7 indicate that there is a 
small discrepancy between the two nonlinear solutions. 
The solution for the more refined model has larger values 
of displacements for some values of the load factor. This 
discrepancy is attributed to the fact that the less refined 
model overestimates the bending stiffness of the shell 
wall. The results in figure 7 also suggest that the 
146 700-DOF model is adequate for conducting nonlin- 
ear analyses for this loading condition. 

Results obtained from nonlinear analyses for a geo- 
metrically perfect shell and for geometrically imperfect 


shells with values for the imperfection-amplitude-to- 
wall-thickness ratio A/t = 0.25, 0.5, and 1 .0 appear in fig- 
ure 8. These results were obtained with the 146 700-DOF 
model and with the imperfection shape in the form of the 
corresponding linear bifurcation buckling mode with a 
negative amplitude. The filled circles in the figure indi- 
cate solutions for the geometrically perfect shell, and the 
unfilled squares, triangles, and circles indicate solutions 
for the geometrically imperfect shells with A/t = 0.25, 
0.5, and 1 .0, respectively. The horizontal dashed line rep- 
resents the linear bifurcation buckling load level. The 
results shown in figure 8 are presented as a plot of the 
load factor p a versus the normal displacement of the shell 
wall at the largest crest of the buckle, defined by the 
coordinates XT = 789.939 in. and 0 = 301.64*, These 
results indicate that the barrel section of the L0 2 tank, 
where the buckles appear, exhibits stable postbuckling 
load-carrying capacity for the geometrically imperfect 
shells and, as a result, is insensitive to initial imperfec- 
tions. Moreover, the shallowness of the barrel panel that 
contains the buckle and the presence of stable postbuck- 
ling load-carrying capacity suggest that the buckling 
behavior is essentially that of the stable-symmetric bifur- 
cation type; that is, the unstable-asymmetric bifurcation 
behavior usually associated with singly curved panels is 
essentially benign. 

The buckle patterns obtained from the linear bifurca- 
tion analysis and the nonlinear analyses are all very simi- 
lar in shape and occur in the same location. The stress 
distribution in the shell that causes the buckle pattern to 
form is shown in figure 9 for a geometrically perfect 
shell and for p a = p b = 1. Contours of the membrane 
meridional and shear stresses on the negative K-axis side 
(0 = 270°) of the tank arc shown in figures 9(a) and 9(b), 
respectively, in units of psi. These contours indicate that 
high meridional compression stresses exist above the 
thrust panel as expected. The contours also indicate that 
buckling occurs in a region of the barrel where there are 
significant shear stresses in the shell that interact with the 
meridional compressive stresses and reduce the buckling 
load. 

Full LH 2 and L0 2 Tanks 

Several different finite-element meshes were used 
for the analysis of the L0 2 tank subjected to the pre- 
launch loading condition with full LH 2 and L0 2 tanks. 
Linear bifurcation buckling analyses were conducted 
first to identify an adequate mesh for modeling the 
behavior of the shell with as few degrees of freedom as 
possible. The passive loads associated with load factor py 
were applied to the STAGS models as a linear prebuck- 
ling stress state (p h = 1), and the active (destabilizing) 
loads associated with load factor p a were used to obtain 
the minimum eigenvalue. 
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The meshes investigated for this prelaunch loading 
condition ranged from 48 990 to 121 500 DOF. The mesh 
identified as adequate for predicting linear bifurcation 
buckling, which appears in figure 10, corresponds to 
99 100 DOF. The reduction in mesh size from 121 500 to 
99 100 DOF was achieved by increasing the mesh in the 
region surrounding the buckle and then eliminating 
unneeded mesh refinement elsewhere in a manner similar 
to that used for the previous loading condition. A 79 950- 
DOF model was investigated that has the same general 
mesh arrangement as that shown in figure 10, but the 
local refinement shown on the right side of the ogive was 
used for both locally refined regions of the ogive. The 
eigenvalues for the 79950- and 99 100-DOF models are 
given by p a ~ 3.94 and p a — 3,78, respectively, where a 
value of 1.0 corresponds to the magnitude of the opera- 
tional loads. Because of the smoothness of the buckling 
mode shown in figure 10 and the 4-percent difference in 
the eigenvalues, the 99 100-DOF model was used to rep- 
resent the linear bifurcation behavior of the L0 2 tank for 
this loading condition. The 79950- and 99 1 00-DOF 
models also were used to obtain nonlinear solutions for a 
geometrically perfect sheik These solutions, which are in 
excellent agreement, indicate that the 99 100-DOF model 
adequately represents the nonlinear behavior of the L0 2 
tank for this loading condition. Thus, all subsequent 
results presented in this section were obtained with the 
99 100-DOF model. 

A short- wavelength buckle in the forward part of the 
aft ogive, shown in figure 10, is essentially a wrinkle in 
the skin on the negative 7-axis side of the tank. The loads 
acting along the shell meridians near the SRB attachment 
point develop the meridional compressive stress result- 
ants shown in figure 1 1 (a) (given in lb/in. for p a = p b -l) 
that caused the buckling mode shown in figure 10. 
Insight into the formation of the buckling mode is 
obtained by noting that Gaussian curvature is a geometric 
measure related to how much the shell membrane stiff- 
ness participates in its bending deformations. The results 
in figure 11 indicate that the meridional compressive 
stress resultants and the circumferential tension stress 
resultants increase in the L0 2 tank from the tip of the for- 
ward ogive to the intertank. The magnitude of the shell 
Gaussian curvature decreases in this part of the shell, 
which causes local buckling to occur at this location. In 
the forward ogive, there are relatively small values of 
hoop tensile stresses to stabilize the shell, but the Gauss- 
ian curvature is high. In the aft ogive, the meridional 
compressive stresses and hoop tensile stresses are larger 
than in the forward ogive, but the Gaussian curvature is 
much smaller. The STAGS results suggest that the higher 
values of meridional compressive stresses and lower val- 
ues of Gaussian curvature in the aft ogive are the domi- 
nant factors that influence the location of the buckling 


mode. These findings also indicate that the shortness of 
the half-wavelength of the buckling mode is a result of 
the high values of the circumferential tensile stress 
resultants shown in figure 11 that are caused by the 
hydrostatic pressure exerted on the shell wall by the L0 2 . 

Results obtained from nonlinear analyses of a geo- 
metrically perfect shell and a geometrically imperfect 
shell with an imperfection-amplitude-to-wall-thickness 
ratio AJt = 0.3 are presented in figures 12 and 13, respec- 
tively. The thickness t in the ratio AJt is the minimum- 
gage wall thickness of the aft ogive and has a value equal 
to 0.1 00 in. The results shown tn figure 13 were obtained 
using an imperfection shape in the form of the corre- 
sponding linear bifurcation buckling mode With a nega- 
tive amplitude. The negative of the linear bifurcation 
buckling mode amplitude obtained was used as the 
imperfection shape because it provided the strongest 
nonlinear interaction with the compressive stresses in the 
shell wall. The results shown in these two figures are for 
nonlinear solutions that were obtained by increasing the 
load factors p a and p b simultaneously to a value of one 
and then holding p b constant while increasing the magni- 
tude of the load factor p a . The normal displacements 
along the length of the aft ogive shell wall are repre- 
sented by the solid lines in the figures for values of the 
load factor p a approximately equal to 3.0, 4.0, and 5.0. 
The linear bifurcation mode is represented by the dashed 
line in the. figures, with normalized amplitudes given by 
the right-hand ordinate of the figures. The bifurcation 
mode in the figures indicates how the imperfection shape 
influences the nonlinear solution. The solid lines shown 
in the figures predict a short-wavelength bending 
response in the aft ogive over the negative 7-axis (0 = 
270°) that is similar in shape to the corresponding linear 
bifurcation buckling mode shape. The overall slope of 
the solid lines (obtained by fitting a straight line to each 
curve) in the figures results from the outward displace- 
ment of the shell wall caused by the increase in pressure 
as p a and p b are simultaneously increased to a value of 
one in the nonlinear analysis. This effect is not repre- 
sented in the linear prebuckling stress state used in a lin- 
ear bifurcation buckling analysis and, as a result, does 
not affect the overall slope of the dashed lines. 

The results presented in figures 12 and 13 predict a 
stable nonlinear response at load levels greater than the 
load predicted by the linear bifurcation buckling analysis 
(p a = 3.78). As the load increases, substantial bending 
gradients (indicated by the waviness of the curves) 
develop and grow in the shell wall, which reduces the 
apparent meridional stiffness of the aft ogive. The non- 
uniformity of the bending gradients is caused by the 
thickness variations in the ogive and the presence of the 
weld lands. These results indicate that a geometrical 
imperfection in the shape of the linear bifurcation 



buckling mode and with a small negative amplitude will 
greatly increase the severity of the bending gradients and 
will cause the growth of the bending gradients to start at 
much lower load levels. 

The reduction in apparent meridional stiffness of the 
aft ogive is shown more explicitly in figure 14. In this 
figure, the intensities of the bending gradients (indicated 
by the magnitude of the normal displacement amplitude) 
at XT = 637.66 in. for the geometrically perfect shell 
(x = 99.3 in. in fig. 12) and at XT = 633.77 for the geo- 
metrically imperfect shell (x = 95.4 in. in fig. 13) are 
given as a function of the load factor p a . These locations 
represent the locations of the largest bending gradients 
shown in figures 12 and 13. The amplitude Aw shown in 
figure 14 is the distance from the maximum value of the 
local shell-wall displacement to the adjacent minimum 
value and represents the intensity of the local bending 
gradient in the response. The filled circles shown in 
the figure correspond to results for a geometrically per- 
fect shell, and the unfilled triangles and squares corre- 
spond to results for geometrically imperfect shells with 
imperfection-ampUtude-to- wall -thickness ratios of A/t = 
0.1 and 0.3, respectively (t = 0.100 in.). The horizontal 
dashed line shown in the figure represents the linear 
bifurcation buckling load level. 

The results shown in figure 14 indicate that the 
amplitude of the greatest local bending gradient grows 
with increasing load and that the amount of growth 
increases substantially with increasing geometric imper- 
fection amplitude. The results predict that the shell can 
support loads greater than the critical buckling load pre- 
dicted by a linear bifurcation buckling analysis. As Aw 
increases, the apparent meridional stiffness decreases, 
and as a result, the positive-valued constant of propor- 
tionality between an increment in load and the corre- 
sponding increment in displacement amplitude 
decreases. This trend is manifested by the reduction in 
slope of the load versus displacement amplitude curves. 
This type of response is similar to the response reported 
by Stevens, Starnes, and Almroth in reference 2 for 
cylindrical shells subjected to combined internal pressure 
and a pure bending moment. The results in reference 2 
indicate that the amplitude of the short-wavelength 
deflection approaches a horizontal tangent as the load 
increases and that the value of the load for the horizontal 
tangent corresponds to a local collapse mode of the cylin- 
der. Mathematically, the horizontal tangent indicates that 
unbounded growth of the displacement occurs for an 
infinitesimal increase in the load. It is expected that the 
curves shown in figure 14 would approach a horizontal 
tangent as Aw increases until a redistribution in load 
occurs within the aft ogive. As a horizontal tangent in a 
load versus displacement amplitude curve is approached. 


the region of the shell containing the bending gradients 
becomes incapable of supporting additional load, and the 
compressive load is redistributed to another portion of 
the aft ogive. If other parts of the ogive cannot support 
the compressive load, the shell will collapse. However, if 
other parts of the ogive can support the compressive 
load, the slope of the load-versus-displacement ampli- 
tude curve will increase as the loading increases. Geo- 
metric imperfections with large amplitudes also could 
cause a similar redistribution in load, and as a result, the 
corresponding slope of the load-versus-displacement 
amplitude curve would increase. Similar curves for 
smaller geometric imperfection amplitudes would 
approach a horizontal tangent. 

The results presented in figures 12 through 14 indi- 
cate that large local bending gradients may occur in the 
shell wall for loads that are much smaller than the local 
collapse load and may cause the TPS to debond from the 
shell wall and fail. The ■. figure 15 

indicate approximate estimates of the local radius of cur- 
vature for the largest bending gradient in the aft ogive 
along the negative F- axis (8 270°), These bending gra- 
dients are located at XT = 637.66 in. for the geometrically 
perfect shell (x = 99.3 in. in fig, 12) and at XT ~ 633.77 
in. for the geometrically imperfect shells (x = 95.4 in. in 
fig. 13). The local radius of curvature p shown in the 
right-hand sketch of the figure was calculated by the 
formula 

2 3/2 

_ [1 +(w') ] 

P “ |w"| 

where w, the local displacement shown in the right- 
hand sketch in figure 15, is approximated by 
w = (Aw/2) $m(2nx/X) where A, is the buckle wave- 
length. The prime marks in the equation denote differen- 
tiation with respect to the local x coordinate. At the crests 
of the wave defined by x = A/4 and 3 A/4, w' = 0, and the 
radius of curvature is given by p = X 2 /(2k 2 Aw). The 
filled circles shown in the figure correspond to results for 
the geometrically perfect shell, and the unfilled triangles 
and squares correspond to results for geometrically 
imperfect shells with imperfection-amplitude-to-wall- 
thickness ratios of A/t = 0.1 and 0.3, respectively. The 
imperfection shape is identical to the linear bifurcation 
buckling mode with a negative amplitude shown in fig- 
ure 13. The results in figure 15 demonstrate that the geo- 
metric imperfection amplitude has a significant influence 
on the local radius of curvature of the shell wall. For 
example, if a given thermal protection system is known 
to debond from the shell wall at a value of p — 100 in., 
the maximum load factor is reduced from a value of 
approximately 4.8 for the geometrically perfect shell to 
3.2 for the geometrically imperfect shell with A/t = 0.3. 
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SWT Full-Scale Structural Test Results 

To assess the accuracy of the STAGS model of the 
SLWT L0 2 tank, finite-element anaylses were per- 
formed on two full-scale structural tests using the same 
modeling approach described previously. The two full- 
scale tests were conducted at the NASA George C. 
Marshall Spaceflight Center on the original SWT during 
the development program of the original Space Shuttle 
ET. Precise measurements of initial geometric imperfec- 
tions in the tank wall were not made for either of these 
test specimens. As a result, these test results can be used 
only as a qualitative means for assessing the accuracy of 
the STAGS model of the SLWT LG 2 tank. 

The SWT L0 2 tank has essentially the same geome- 
try as the SLWT L0 2 tank, but is made of 2219 alumi- 
num alloy. The primary difference between the two L0 2 
tanks is that the skins of the SWT are thicker than those 
of the SLWT, with the thicknesses much more uniformly 
distributed over the SWT shell. Thus, the modeling 
approach described for the SLWT was used to model the 
two full-scale SWT test articles. The SWT model was 
generated by modifying the STAGS user- written subrou- 
tine WALL for the SLWT to account for the SWT thick- 
nesses. The differences between the SWT intertank and 
the SLWT intertank are negligible for the purpose of 
transferring loads from the intertank to the L0 2 tank. The 
two full-scale SWT tests described subsequently are 
referred to herein as the structural test article (STA) and 
the ground vibration test article (GVTA). 

STA Results 

The STA consisted of a SWT L0 2 tank and a SWT 
intertank mounted vertically to a LH 2 tank load simulator 
and two rigid vertical posts at the SRB attachment points. 
The LH 2 tank load simulator was modeled by a self- 
equilibrated line load applied to the bottom of the inter- 
tank, as indicated in figure 16. A uniformly distributed 
circumferential line load of 1394 kips was also applied to 
the tank at XT — 852.8 in. For the test, these two loads 
were applied and then the tank was filled with room tem- 
perature water while an ullage pressure was maintained 
in the tank. After filling the tank to XT = 455 in., a depth 
of approximately 42 ft, the ullage pressure was slowly 
reduced. When the ullage pressure reached 0.57 psig, the 
tank unexpectedly buckled in the forward ogive between 
XT = 455 in. and XT = 475 in. and between 0 = 253° and 
277° (negative F-axis side of the tank). The SRB inter- 
face forces that are reacted at the two vertical posts had 
magnitudes equal to 1295 kips (fig. 16). 

To simulate the test loading conditions in a practical 
manner, all loads shown in figure 16, except for the SRB 
beam loads, were controlled by the load factor p a . A 
value of p a = 1 corresponds to values of the loads at 


which the actual test article buckled. The loads at the 
ends of the SRB beam were computed as reactions and 
found to be in excellent agreement with the correspond- 
ing SRB beam loads shown in figure 1 6. 

Finite-element meshes were constructed for the STA 
and a limited convergence study was performed that fol- 
lowed an approach similar to that described for the 
SLWT analyses. The final mesh used to analyze the STA 
has 159993 DOF and is highly refined in the forward 
ogive on the negative F-axis side of the tank. The linear 
bifurcation buckling mode obtained for the STA with this 
mesh is shown in figure 17. This buckling mode is a 
short- wavelength buckle similar to the one obtained for 
the SLWT prelaunch loading condition with full LH 2 and 
L0 2 tanks. The location of the buckling mode shown in 
this figure is the same as the location observed during the 
test. The eigenvalue is given by p si = 1.14. 

Next, a series of nonlinear analyses were conducted 
using the 159993-DOF mesh for values of the 
imperfection-amplitude-to-wall-thickness ratio A/t = 0, 
0.1, 0.25, 0.5, and 1.0. The thickness t in the ratio A/t is 
the minimum-gage wall thickness of the forward ogive 
and has a value equal to 0.080 in. For each of these cases, 
the geometric imperfection shape was input in the form 
of the linear bifurcation buckling mode shown in fig- 
ure 17, with a negative amplitude to obtain the strongest 
interaction between the membrane compressive stress 
and shell-wall deformations. 

Figure 18 shows the results of the nonlinear analyses 
of the STA. The maximum normal displacement that 
occurs at the crest of the buckle pattern is shown in the 
figure as a function of the load factor p a . The buckle crest 
is located at XT = 457.6 in. and XT = 466.6 in. for the 
geometrically perfect and imperfect shells, respectively, 
and at 0 = 267.2°. The filled circles and the unfilled cir- 
cles, diamonds, squares, and triangles correspond to 
results for A/t = 0, 0,1, 0,25, 0.5, and 1.0, respectively. 
The results for all values of A/t indicate a monotonic 
increase in load with increasing normal displacement. 
However, a maximum load was reached at which numer- 
ical difficulties were encountered in the nonlinear solu- 
tions. After several unsuccessful tries to increase the load 
above these values, it was concluded that the last data 
point on each curve corresponds to a limit point of the 
shell response. At each limit point, the shell buckled into 
a mode similar in shape to the buckling mode shown in 
figure 17 and at the same location. The meridional and 
circumferential stress resultant distributions are similar 
to those presented in figure 11 for the prelaunch loading 
condition with full LH 2 and L0 2 tanks. 

The values of the limit points (filled circles) are 
shown in figure 19 as a function of geometric imper- 
fection amplitude and load factor, and indicate the 
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imperfection sensitivity of the STA. The limit points 
obtained for the STA span a broad load range bounded 
by P a ~ 1*18 for a geometrically perfect shell and 
p a ~ 0.53 for a geometrically imperfect shell with 
A/t - 1.0 (A - 0.080 in.). This load range corresponds to a 
55-percent reduction in load-carrying capacity of a geo- 
metrically perfect shell. The dashed line shown in the 
figure has a value of A/t = 0.064 for p a = 1 (which corre- 
sponds to the buckling load of the test). This result sug- 
gests that the STAGS modeling approach provides a 
reasonably accurate indication of the SWT behavior. 
Because of the similar geometric character of the SWT 
and the SLWT, the results also suggest that the SLWT 
model should provide a reasonable representation of the 
SLWT nonlinear shell response. 

GVTA Results 

The GVTA consisted of a SWT mounted on two 
SRBs arid an orbiter attached to the SWT. The SWT con- 
sisted of a L0 2 tank, a LH 2 tank, and an intertank. In this 
configuration, the SWT is inclined at an angle of approx- 
imately 10° in the XT-Z plane because of the eccentric 
weight of the orbiter. The loads acting on the L0 2 tank 
and intertank during the test and the inclination angle a 
are shown in figure 20. These loads consist of two SRB 
interface force components, the LH 2 tank interface force 
and moment, a uniformly distributed circumferential line 
load of 20.86 kips applied at XT — 852.8 in., and a hydro- 
static water pressure distribution that corresponds to the 
tank fill level of XT = 645 in. (a depth of approximately 
26.5 ft). No ullage pressure was present inside the L0 2 
tank during the test. The hydrostatic pressure distribution 
for the GVTA was defined in the STAGS model with 
user-written subroutine UPRESS in terms of the local 
axial coordinate x shown in figure 20(b). The pressure 
distribution is given by p(x, 0) = 0 for values of 
x<Xf — r(x) tan 01 cos0, where Xf is the local 
coordinate of the fill level defined by the positive 
numerical difference between stations XT - 645 in. and 
XT - 371 in., as shown in figure 20(b). The symbol 0 is 
the cylindrical coordinate defined in figure 2, and r(x) is 
the polar radius of the shell reference surface (fig. 20(b)) 
that was calculated from the differential geometry of the 
L0 2 tank components. For the remaining values of x, the 
pressure is given by 

p(x, 0) = y[(x-Xy) cosa + r(x) sina cos0] 

where y represents the specific weight of water at room 
temperature. 

The original test plan for the GVTA was to fill the 
tank with water and then to perform a ground vibration 
test. However, when the water level reached XT = 645 
in., the tank unexpectedly buckled in the forward ogive 


between approximately XT = 437 in. and XT = 503 in. 
and between 0 = 247° and 281° (negative F-axis side of 
the tank). 

The test loading conditions for the GVTA were sim- 
ulated in the manner described herein for the STA; that 
is, all loads shown in figure 20, except for the SRB inter- 
face forces, were assigned to the load factor p a such that 
a value of p a = 1 corresponds to buckling of the test arti- 
cle. The loads at the ends of the SRB beam were com- 
puted as reactions and found to be in excellent agreement 
with the SRB interface forces shown in figure 20(a). 

Several finite-element meshes were also investigated 
for the GVTA following the same approach used for the 
STA. The final mesh selected to analyze the GVTA has 
252 300 DOF and is highly refined in the forward ogive 
on the negative F-axis side of the tank. The linear bifur- 
cation buckling mode obtained for the GVTA with this 
mesh is shown in figure 21. This buckling mode is also a 
short- wavelength buckling mode similar to the one 
obtained for the SLWT prelaunch loading condition with 
full LFI 2 and L0 2 tanks and for the STA. The location of 
the buckling mode shown in this figure is essenfi&Hy the 
same as the location observed during the test. The eigen- 
value is given by p a = 2.41 . 

Next, the 252 300-DGF mesh was used to conduct 
nonlinear analyses for values of the imperfection- 
amplitudc-to-wall-thickness ratio A/t = 0, 0.125, 0.25, 
0.5, 0.75, and 1.0. The minimum-gage wall thickness t 
for the forward ogive has a value equal to 0.080 in. For 
each of these cases, the geometric imperfection shape 
was in the form of the linear bifurcation buckling mode 
shown in figure 21 with a negative amplitude. 

The results of the nonlinear analyses of the GVTA 
are shown in figure 22. The maximum normal displace- 
ment that occurs at the crest of the buckle pattern is 
located at XT = 466.7 in. and at 0 = 267.2°, and is shown 
in the figure as a function of the load factor p a . The filled 
squares, triangles, and circles and the unfilled circles, 
squares, and triangles correspond to results for A/t = 0, 
0.125, 0.25, 0.5, 0.75, and 1.0, respectively. Similar to 
the STA results, the GVTA results indicate a monotonic 
increase in load with increasing normal displacement that 
terminates at a limit point for all values of A/t. At each 
limit point, the shell buckles into a mode similar in shape 
to the mode shown in figure 21, and at the same location. 
The values of the limit points (filled circles) shown in 
figure 23 as a function of the geometric imperfection 
amplitude and load factor indicate the imperfection sen- 
sitivity of the GVTA. The limit points obtained for the 
GVTA span a broad load range bounded by p a = 2.46 for 
a geometrically perfect shell and p a - 0,92 for a geomet- 
rically imperfect shell with A/t = 1 . This load range cor- 
responds to a 62-percent reduction in the load-carrying 
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capacity of a geometrically perfect shell. The load reduc- 
tion for the GVTA is approximately 7 percent greater 
than that for the STA. The dashed line shown in figure 23 
has a value of AJt = 0.625 (A = 0.050 in.) for p a = 1 
(which corresponds to the buckling load of the test). This 
result also suggests that the STAGS modeling approach 
provides a reasonably accurate indication of the SWT 
behavior and the SLWT nonlinear shell response. 

Concluding Remarks 

The results of buckling and nonlinear analyses of the 
Space Shuttle superlightweight tank (SLWT) liquid- 
oxygen (L0 2 ) tank have been presented. An overview of 
the L0 2 tank and intertank structures and the loading 
conditions for two important prelaunch loading condi- 
tions have been described. In addition, the analysis code 
used in the study has been described and the finite- 
element modeling approach and details presented. The 
analytical method used in the study to simulate the load- 
ing conditions associated with prelaunch fueling of the 
Space Shuttle has been discussed. 

Results have been presented for the superlightweight 
L0 2 tank subjected to two prelaunch loading conditions 
and for two full-scale structural tests that were conducted 
during the development program of the original standard- 
weight tank (SWT). These results illustrate three dis- 
tinctly different types of nonlinear response for thin- 
walled shells subjected to combined mechanical and 
thermal loads that may be encountered in the design of 
other liquid-fuel launch vehicles. Predicting the response 
of these shells generally requires large-scale, high-fidel- 
ity finite-element models to represent the response accu- 
rately. For the first SLWT prelaunch loading condition, 
the liquid-hydrogen (LH 2 ) tank is full and the L0 2 tank 
is empty. The analytical results predict that the nonlinear 
response is characterized by a buckling response that is 
insensitive to initial geometric imperfections. For this 
loading condition, the barrel section of the L0 2 tank is 
predicted to buckle at loads that are more than twice the 
operational loads. For the second SLWT prelaunch load- 
ing condition, the LH 2 and L0 2 tanks are full. The 
nonlinear response for this loading condition is charac- 


terized by a short-wavelength bending gradient that 
grows in amplitude in a stable manner with increasing 
load. For this loading condition, local bending gradients 
appear in the aft ogive of the L0 2 tank that do not lead to 
a general instability mode but may cause failure of the 
thermal protection system for load levels in excess of 
approximately twice the operational load level. More- 
over, the results predict that the severity of the local gra- 
dients is significantly affected by localized initial 
geometric imperfections. 

For the two full-scale structural tests of the SWT, the 
nonlinear responses exhibit local buckling of a doubly 
curved shell segment of the L0 2 tank forward ogive that 
is characterized by a limit-point behavior. The magnitude 
of the load level corresponding to the limit point has 
been shown to be very sensitive to local initial geometric 
imperfections in the L0 2 tank. Specifically, load reduc- 
tions of about 55 to 62 percent of the buckling load of a 
geometrically perfect shell are predicted for a geometric 
imperfection shape in the form of the linear bifurcation 
buckling mode and with a one-wall-thickness imperfec- 
tion amplitude. The buckling loads obtained from both 
tests correspond to geometric imperfection amplitudes 
that are less than one minimum-gage wall thickness. For 
both tests, the analytical results suggest that the finite- 
element modeling approach used in the present study 
represents the nonlinear behavior of the superlightweight 
L0 2 tank very well. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
September 18, 1996 
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Figure 4. Loading characteristics. 



Figure 5, Axisymmetric temperature profile for prelaunch loading condition with full LH 2 and L0 2 tanks. XT' values in 
inches. 
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Figure 6. Linear bifurcation buckling mode for prelaunch loading condition with full LH 2 tank and empty L0 9 tank 
146 700 DOF; XT values in inches. 


104 600-DOF model 




Figure 7. Convergence of nonlinear solutions for prelaunch loading condition with full LH 2 and empty L0 2 tanks. 
Imperfection-amplitude-to- wall-thickness ratio AJt = 0.25; XT value in inches. 
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Figure 8. Effect of imperfection amplitude on nonlinear solutions for prelaunch loading condition with full LH 2 and 
empty L0 2 tanks. XT value in inches. 
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Figure 9. Membrane stress distributions in geometrically perfect L0 2 tank and intertank for prelaunch loading condition 
with full LH 2 and empty L0 2 tanks. p a =p b = 1; stress values in psi. 
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Figure 12. Bending gradients in aft ogive of geometrically perfect L0 2 tank shell for prelaunch loading condition with 
full LH 2 and L0 2 tanks for & - 270°. XT value in inches. 
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Figure 13. Bending gradients in aft ogive of geometrically imperfect L0 2 tank shell for prelaunch loading condition 
with full LH 2 and L0 2 tanks for 0 = 270°. Imperfection-amplitude-to-wall-thickness ratio A/t = 0.3; XT value in 
inches. 
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Figure 14. Effect of imperfection amplitude on nonlinear solutions for prelaunch loading condition with full LH 2 and 
L0 2 tanks for 0 = 270°. XT values in inches. 
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Figure 15. Effect of imperfection amplitude on local radius of curvature of largest bending gradient for prelaunch load- 
ing condition with full LH 2 and LQ 2 tanks for 0 = 270°. XT values in inches. 
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Figure 16. Loads for SWT STA at buckling. XT values in inches. 



Figure 17. Linear bifurcation buckling mode for STA, 159 993 DOF. XT values in inches. 
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Figure 18. Effect of imperfection amplitude on collapse load for STA. XT values in inches. 
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Figure 19. Predicted imperfection sensitivity of STA. 







(b) Side view. 


Figure 20. Loads for SWT GVTA at buckling. XT values in inches. 





Figure 21 . Linear bifurcation buckling mode for GVTA. 252 300 DOF; XT values in inches. 
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Figure 22. Effect of imperfection amplitude on buckling load for GVTA. XT values in inches. 
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